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Eigen modes for the problem of anomalous light transmission 
through subwavelength holes 
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M. Gorkunov 
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We show that the wide-spread concept of optical eigen modes in lossless waveguide structures, 
which assumes the separation on propagating and evanescent modes, fails in the case of metal- 
dielectric structures, including photonic crystals. In addition to these modes, there is a sequence 
of new eigen-states with complex values of the propagation constant and non-vanishing circulating 
energy flow. The whole eigen-problem ceases to be hermitian because of changing sign of the 
optical dielectric constant. The new anomalous modes are shown to be of prime importance for the 
description of the anomalous light transmission through subwavelength holes. 
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The concept of eigenmodes lies at the center of any 
physical problem for optical light-transmitting systems 
including waveguides and photonic crystals. This con- 
cept seems to be clear and the general properties of the 
eigen-values and eigen-functions seem to be fully pre- 
sented in numerous textbooks, see, e.g., Q, 0,0, Of. 

For photonic crystals (PC), determination of the band 
structures - the eigen-frequencies and eigen-functions 
versus the wavevector - is the most common. The eigen- 
frequencies are proven to be real for any lossless PC, i.e., 
for any periodic distribution of real dielectric optical per- 
mittivity e(r). A close analogy with the electronic band 
states is well established 

Another formulation of the eigen-mode problem is typ- 
ical of waveguiding systems which are uniform along the 
propagation coordinate z: To find the admitted values 
of the propagation constant (3 entering the propagation 
factor exp(i/3z) and the corresponding transverse (in x 
and y) distributions of the optical light fields. 

It is known for the lossless dielectric systems (e > 0) 
that the admitted values of (i 2 are real - positive and neg- 
ative 0, 0] • The corresponding eigenmodes are propa- 
gating and evanescent, respectively, and a close analogy 
with quantum mechanics holds true. The same is valid 
for the ideal- metal waveguiding systems, e — > oo. 

A permanently growing interest to metal-based sys- 
tems is usually attributed to the excitation of surface 
plasmons [l], |5| which allow the light constrain on a sub- 
wavelength scale. The necessary conditions, e' — Ree < 
— 1, e" — Ime <C 1, are fulfilled for many metals. The 
limit of lossless metal (e' < 0, e" = 0) is as actual as 
that of lossless dielectric. Recent discovery of the ex- 
traordinarily light transmission through subwavelength 
holes in metal slabs ||g, |6| has strongly enhanced the cur- 
rent interest. However, the main ingredients and para- 
metric dependences of this fundamental phenomenon re- 
main poorly understood; its description relies on numer- 
ical simulations and simplified models 0, El . 

The present notion of eigenmodes for lossless (e" = 0) 
metal-dielectric waveguide structures (ID and 2D PCs, 
arrays of holes and single holes in metals, etc.) is sur- 
prisingly scarce. It is widely believed that the values of 



(3 2 remain real, i.e. the separation into the propagat- 
ing and evanescent modes and the analogy with quan- 
tum mechanics, still hold true. Using particular mod- 
els for e(x,y) and Fourier expansions, many authors cal- 
culate numerically only the real values /3 2 (k x , k v ), see, 
e.g., 0, E| and references therein. To the best of our 
knowledge, there are only few publications where the au- 
thors mention, on the basis of particular calculations, 
that some of values of 1 cease to be real 0, . 

The aim of this letter is {%) to show that the gener- 
ally accepted concept of eigenmodes is insufficient for 
the case of lossless metal-dielectric waveguide structures, 
(ii) to analyze the general properties of new eigenmodes 
which are neither propagating nor evanescent, and (Hi) 
to demonstrate that these new modes are crucial for the 
description of the anomalous light transmission through 
subwavelength holes. We will also pay attention to the 
propagating modes in the subwavelength case. 

To catch the essence of the eigen-mode problem, we 
consider the case of a ID PC with e being a real periodic 
function of the coordinate x. Restricting ourself to the 
TH polarization, we have from Maxwell equations for the 
only nonzero component of the magnetic field H = H y : 



LH = f3 2 H , 



L- d 1 d 

dx e dx 



ek 



(1) 



where fco = 2tt/\. This sets an eigen-mode problem 
for the operator L with 1 being the eigen-value. The 
nonzero components of the light electric field are ex- 
pressed by H: E x = (3H/ek , E z = (i/ek ) dH/dx. 

The conclusion about reality of (3 2 and analogy with 
quantum mechanics is generally based on Hermitian char- 
acter of the differential operator L, i.e., on the prop- 
erty (H 2 \LH l ) = (H 1 \LH 2 )*, where (..|„) stands for the 
scalar product. The only possibility to get this property 
is to define the scalar product of two complex functions 
Hi and H 2 as 



^(.t) H;(x)H 1 (x) dx . 



(2) 



In the dielectric case, s(x) > 0, this definition meets all 
axioms of the scalar product [14( . For sign-changing func- 



2 



tions e(x), it contradicts, however, to one of the axioms - 
the norm (H\H) ceases to be positively defined. There- 
fore, the conclusions about Hermitian character of L and 
about reality of 1 cannot be made. 

Apart from the above negative observation, some pos- 
itive assertions can be made: 

- If an eigen-function possesses a positive norm, the cor- 
responding value of 1 is real. Two eigen-functions Hi 
and H 2 corresponding to different real eigenvalues 01 and 
02 are orthogonal, (i/^l-Hi) = 0. If at least one of /3i,2 is 
complex, it has to be replaced by (H 2 \H\) = 0. 

- In the case of TE polarization (nonzero E yi H x , and 
H z components) the above peculiarities are absent - the 
admitted values 1 are real and an analogy with quan- 
tum mechanics holds true. 

- In the 2D case, Hermitian character of the eigen-mode 
problem and reality of 1 can be proven in exceptional 
cases of waves with E z = (TE modes). In the general 
case (TM and hybrid modes) complex values of /3 2 , i.e., 
anomalous modes, are present, see also below. 

Now we turn to particular examples. Consider first 
a ID PC consisting of an alternating sequence of lay- 
ers with permittivities £1,2 and thicknesses x\$\ the pe- 
riod is x = x i + x 2- Solution of Eq. (1) for an eigen- 
mode has the form H(x) — e lnx h(x), where K is the 
Bloch wavevector ranging from — w/xq to n/xo and h(x) 
is an a;o-periodic function. Using continuity of H(x) and 
E z (x) oc e~ 1 dH/dx at the interfaces and periodicity of 
h(x) and dh/dx, we obtain the followingdispersion equa- 
tion for 1 (it generalizes Eq. (13) of [l2|'): 



1 / Pi£2 Vi£\ 



sin(pixi) sin(p 2 x 2 ) 



2 \p 2 £i P1E2/ 

= cos(pia;i) cos(p 2 x 2 ) — cos(«;a;o) 



(3) 



where pi 2 = (£1,2^0 ~~ 2 ) 1 ^ 2 - The trigonometric func- 
tions possess generally complex arguments. If £i ;2 are 
real, Eq. (3) is real as well. The roots j3 = iiy/ffi de- 
scribe the wave propagation in the ±z directions. 

The circles in Fig. la show the values of the propa- 
gation constant on the complex 0,0" plane in the loss- 
less case for k — (at the center of the Brillouin zone), 
x i = x /b = A/4, £1 = 1, and £ 2 = —9.6. The last pa- 
rameter corresponds to silver at A w 0.5 /im [l^. We 
have a single real root, 0/ko « 1.2, and two sequences of 
roots with 0' 7^ 0. One sequence, with pure imaginary 
values of refers to the known evanescent modes. An- 
other sequence, with 0^0 and 0' > \/\e 2 \ (complex 
values of 2 ) corresponds to anomalous modes forbidden 
in the dielectric case. They combine features of prop- 
agating and evanescent modes. The symmetric pairs of 
anomalous roots correspond to mutually conjugate values 
of/3 2 . 

At k = 0, the modes are either even, h(x) = h(—x), 
or odd, h(x) = —h(—x); the "even" and "odd" roots are 
shown by filled and open circles in Fig. la. The single real 
root corresponds to the even propagating mode. About 
30% of the modes are anomalous. 





c 

• 


• 








( 
1 

c 


5 
1 

) 








1 

( 
« 


C 

» 








• j 

! 










x, = M 

x 2 = ^ 




C 

(a) 


> 

























i Ht- 


1 


W 






-» 




O 






O 






• 


• 






( 


C 

0) 



























, — 
■ 


■ 










■ 
■ 






m 






■ 
■ 








■ 

1 


■ 

1 








| Ft = X/2 


« 


1 


1 





















-0.2 0.0 0.2 1.2 

P'/k„ 



-0.1 0.0 0.1 1.2 



-0.02 0.00 0.02 0.80 0.82 



FIG. 1: Values of j3 for £2 = —9.6 and three structures: 
(a) ID PC with X2 = 4a; 1 (the gray line shows positions of 
evanescent roots for e'2 = 0.3); (b) a single slit with xi = A/4; 
(c) a single cylindric hole with R = A/2. The vertical solid 
lines in (b) and (c) show the continuous spectrum. 



With weak losses taken into account, the roots experi- 
ence small displacements. The strongest of them (shifts 
to the right) occur for the evanescent roots. The gray 
line in Fig. la shows positions of these roots for e 2 ' = 0.3, 
which models silver at w 0.5 /xm. The other displace- 
ments do not exceed the circle size. 

Whereas the Pointing vector P = (P X ,P Z ) is zero for 
the evanescent modes, it is nonzero for the anomalous 
modes. Each anomalous mode is responsible thus for cir- 
culation of the light energy. Applying the orthogonality 
relation (i?||ifi) = to a conjugate pair of anomalous 
modes, we come to the following relation for an anoma- 
lous eigen-function h a (x): 



-\x) \h a (x)\ 2 dx = 



(4) 



where the integration occurs over a period of Xq. Since 
P z cx 0e- 1 (x)\h(x)\ 2 exp(-20'z), it shows that the x- 
averaged energy flux is zero for the anomalous modes - 
the inflow in air is compensated by the outflow in metal. 
Eq. J2J can be interpreted also as a zero norm for the 
eigen- vector h a (x). The eigen- vectors h(x) for the prop- 
agating and evanescent modes possess positive norms. 
Fig. 2 shows the energy circulation within a half-period 
Xq/2 for the first anomalous mode with > 0. For the 
mode with < the arrows have to be inverted. 

What happens with the above modes when changing 
2:1,2? With decreasing slit width x\, the propagating 
mode survives (0 grows as x^ ), the evanescent mode 
with 0" < v|£2| disappears, the other evanescent modes 
experience minor changes, and the vertical distance be- 
tween the anomalous modes increases as 2tt/x\. 

Increasing wall parameter x 2 does not affect strongly 
the propagating and anomalous modes in Fig. la. How- 
ever, the density of evanescent modes in the region 0' > 
•v/|£ 2 | increases oc x 2 - The limit x 2 — » 00 corresponds 
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FIG. 2: Streamlines of the Pointing vector for the first anoma- 
lous symmetric mode with /3' > 0. 



to the single-slit case. It is represented by Fig. lb. We 
have single localized (in \x\) propagating and evanescent 
modes, a sequence of localized anomalous modes, and a 
continuous spectrum of non-localized evanescent modes 
with j3' = and j3" > \f\s2\- This also follows from a 
direct description of the single-slit case. The anomalous 
modes represent thus properties of a single-hole. 

For koX\ <C 1, X2 — > oo (very narrow single slit) we 
have for the single propagating mode and e 2 ' <C |e 2 |: 



Xi 



In 



1 



1 



2% 4 



1 



(5) 



The effective refractive index (3' /kg tends to infinity for 
xi — > 0. In the ideal-metal limit ( | er 2 1 — * oo) we have 
(3 = 0. For x±.2 — * oo the propagating mode transforms 
to the usual surface plasmon. 

The eigen-functions for the above periodic and single- 
slit cases (localized and delocalized) can be constructed 
readily from the exponential functions exp(zp 12 a;) with 
proper values of pi i2 (j3 2 ). 

The Bloch dependences (3'(k) and (3"(k) for the prop- 
agating and evanescent modes, which are present in the 
PC literature, are insufficient in the metal-dielectric case 
where the anomalous modes are present. Moreover, the 
anomalous modes cause strong interaction of different 
branches and peculiarities of Bloch diagrams. Fig. 3 
shows examples of the ^-dependences for our ID PC 
structure. The value of \/3'\ for the first pair of anoma- 
lous modes (see Fig. la) decreases with growing k and 
turns to zero for 0.66 < kxq/t: < 0.71, i.e., the anoma- 
lous modes become evanescent within this interval. For 
kxq/tt > 0.71 a pair of anomalous modes appears again. 
This behavior is correlated with bifurcations of /?"(«) 
for the anomalous modes and the nearest evanescent 
mode. Behavior of the second pair of anomalous modes 
is different. The confluence of (3'(k) for these modes at 
kxq/tt ~ 0.9 is accompanied by the split of /?"(«). The 



FIG. 3: Dependences /?"(«) and /3'(k) for the periodic struc- 
ture with xo = 5xi — 5A/4, and £2 = —9.6. Solid lines are 
plotted for the first pair of anomalous modes and the nearest 
evanescent mode, the dotted lines correspond to the second 
anomalous pair, and the gray line is for the propagating mode. 



dependence @(k) for the propagating mode and most of 
evanescent modes is fairly weak. 

The situation described is also inherent in 2D struc- 
tures. Consider for simplicity TM (H z = 0) modes for a 
cylindric hole of radius R in a lossless metal. The corre- 
sponding dispersion equation reads |l6j | 



P2 Jl #0 (1) - £2Pl Jo H[ 1} = , 



(0) 



where Jo,i = Jo,i(pi-R) and = Hq 1 1(p 2 R) are the 

Bessel and Hankel functions. Fig. lc shows the spectrum 
of f3 for £2 = —9.6 and R = A/2. One propagating mode, 
two evanescent modes, a sequence of localized anomalous 
modes, and a continuous spectrum are present. 

The TM propagating mode exists here only for R > 
Rmin{£2) with Rmin decreasing from ~ 0.334A to ~ 
0.257A when £2 changes from —10 to —1. For the HEn 
propagating mode, which provides the minimum cut-off, 
R m in decreases from ~ 0.22A to ~ 0.087A in this range 
of £2- For the ideal metal R m i n /\ ~ 0.3. As follows from 
Eq. |JBJ, in the range — 1 < £2 < a propagating mode 
exists even for R — > 0. 

To show the importance and applications of our ap- 
proach, we consider the following fundamental problem: 
An x-polarized plane wave of a unit amplitude is incident 
normally onto the above ID metal-dielectric PC occupy- 
ing the semi-space z > 0. It is necessary to determine 
the transmission properties at the interface z = 0. 

We use the following explicit mode expansion for H: 

§ a v h v {x) e*^ z ; z > 

"=< ^ (7) 

e ik Q z + g bn e inKx + q n z . z<Qj 
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FIG. 4: Dependence T(X/xo) for £2 = —9-6, 2:2/2:1 = 4, and 
iV = 48; the solid and dashed lines are plotted for e 2 ' = and 
0.3. The inset shows T(\/xo) for e'2 = and five values of 
X2/X1 in the vicinity of the main peak. 



where v numerates the even eigen-modes at k = 0, 
K = 2tt/x , q n = (n 2 K 2 - k 2 ) 1 ' 2 for n 2 K 2 > k 2 , 
q n = -i{k 2 - n 2 K 2 ) x l 2 for k\ > n 2 K 2 , and b n and 
a v are the amplitudes to be found. Real and imaginary 
values of q n correspond to evanescent and propagating 
waves in air. With increasing period, new diffraction or- 
ders appear at xq = n\ in reflected light; only a single 
propagating wave (with n — 0) is allowed for x < A. 

Employing the boundary conditions at the interface 
z = and the Fourier transformation, we come to the 
set of linear algebraic equations for determination of a„: 

00 

J2 au[K - * &«n 1 (h"/e) n ] = 2S n0 , (8) 

v=0 

where h v n and [h v je) n are the Fourier components of the 
periodic functions h v {x) and h u (x)/e(x). These compo- 
nents were expressed explicitly through (3 V . 

A numerical routine with truncation of this set at 
v-max = n max = N 3> 1 was used to calculate a v and de- 
termine the reflection and transmission properties. The 
results of calculations converge perfectly well with in- 
creasing N when all eigen-modes are taken into account. 
However, the calculation results diverge when the anoma- 
lous modes are ignored. This means (i) that the set of 
the eigen-functions is incomplete without the anomalous 



modes and (ii) that these modes are directly involved in 
the energy transfer from the incident plane wave to the 
propagating mode(s) inside the slits. 

Consider the key results for the transmission coeffi- 
cient T defined as the energy-flax fraction transmitted 
into the propagating mode at z = 0. This definition is 
useful even for e" =/= when the propagating mode de- 
cays inside the metal as exp(— 2/3"z). Fig. 4 shows the 
dependence T(A) for the lossless and lossy cases. The 
losses influence T very weakly; their main effect occurs 
during propagation. The function T(X) possesses narrow 
peaks at A = xo/n; the main peak corresponds to n = 1. 
The geometric nature of the resonances is evident and 
it is not overshadowed by the Fabry-Perot interference 
typical of the slab case Q • The maximum value of the 
transmission coefficient, T max ~ 0.47 (which corresponds 
to x\ = A/5) is pretty high. It is decreasing slowly with 
increasing |e 2 |- For A > xq, the difference 1 — T is the 
reflection coefficient. 

Our procedure allows to resolve the fine structure of 
the main peak and investigate the dependence T max on 
the ratio x^/xi, see the inset in Fig. 4. The structure 
is different for x 2 /xi < |ea| and x 2 jx\ > |ea | . In the 
first case, the peak spike is sharp and centered exactly 
at A = Xq. In the second case, the right wing of the 
peak is broadened and the very maximum is slightly 
shifted to the right. Surprisingly, T max {xi/x2) is not 
a monotonously increasing function; the optimum value 
of x\jxi is ~ 1/6 where T max ~ 0.51. Note that the 
effective permittivity of our medium in the limit A Xq 
is (e -1 ) -1 = MW^il^l ~x 2 ), see, e.g., H3; it turns 
to infinity and changes sign at x<2,jx\ = \e<^. Validity of 
this limiting case was controlled numerically. 

In conclusion, the eigen-mode problem for metal- 
dielectric structures - photonic crystals, single holes in 
a metal - is essentially different from that typical of di- 
electric structures and quantum mechanics. It is not 
hcrmitian and the eigen-values for the localized states 
are generally complex in the lossless case. The anoma- 
lous transmission through subwavelength holes is sub- 
stantially determined by the transmission coefficient for 
a single interface. This characteristic possesses a num- 
ber of fundamental features. The opposite faces of the 
metal-dielectric structure arc coupled via a weakly de- 
caying propagating mode. 
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